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Q I The thermal and quark/baryon chemical potential dependences of quark condensate and masses 

O ■ of IT- and (T-mesons are studied in the instanton model of the QCD vacuum in precritical region. 

^ The impact of phonon-like excitations of instanton liquid on the characteristics of a-meson in 
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such an environment is also examined. 



Nowadays with the common belief that the hadron world and its symmetries are driven by QCD 
our understanding of the lightest mesons origin, at least from the hadron side, is tightly related to the 
t> I spontaneous breakdown of chiral symmetry inherent to the corresponding Lagrangian. Then the pions 
appear as the massless Goldstone bosons while the scalars acquire a non-zero vacuum expectation 
values reflecting the sophisticated structure of QCD vacuum. At the level of QCD fundamental flelds 
^ . we have phenomenologically learnt that this symmetry breaking is dominated by the non-zero values 
[ of quark condensates [1] and the scalar meson attributes result from the dynamical generation of light 
O [ quark masses. However, at rather high temperature T, as the lattice QCD studies have taught [2], the 
^ ■ phase transition restoring the chiral symmetry occurs (and a quark chiral condensate is just its order 
^ parameter) together with the transition deconflning the quarks and gluons and another transition of 
^ colour superconductivity (characterized by a diquark condensate) at high quark/baryonic chemical 
r-| ! potential /i. Clearly, the last phase transition is scarcely feasible in the terrestrial conditions and is 
j> I relevant rather in the astrophysics observations of compact stars. As to two others they are intensively 
] explored in ultrarelativistic heavy ion collisions [3]. 

■ The properties of various phases of hot and dense matter are dependent on the interplay 
between the competing processes of quark-quark and quark-antiquark pairing (influenced by the 
QCD vacuum) and, as expected, could be clarifled by studying the QCD phase diagram on (/i— 
T)-plane in the lattice simulations (albeit a very idea is relevant for the infinite matter and in 
thermodynamic equilibrium). The first results of lattice phase diagram studies were quite schematic 
but, at the same time, indicative enough because they generated very important question about the 
structure of phases around the deconfinement critical temperature Tc, i.e. at \T — Tc\ -C Tc. Today as 
become especially clear when it was found that J/ ip charmonium survives at rather high temperatures 
beyond Tc [4j. In fact, this result has been confirmed by the RHIC experimetal observations and has 
led to the concept of strongly coupled quark-gluon plasma [5j revealing the existence of enormous 
number of bound states at T > in deconfined phase. 

Phenomenological analysis of strongly interacting matter under the extreme conditions is 
complicated by the comparatively low credibility of dynamical hadron models (linear and nonlinear 
sigma-models, the Nambu-Jona-Lasinio model and its modificatios etc.) on the (yU - T)-plane. Their 
predictions are not likely to be reliable at approaching the critical line of phase transition restorating 
the chiral symmetry. At small values of and T the parameters of those models are tuned in 
such a way to reproduce some experimental data and to obey some general constraints rooted in 
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QCD. However, extending the models to the phase transition regions leads sometimes to conflicting 
conclusions (see, for example, the recent discussion in ^) At the same time the analysis of precritical 
regions could be entirely informative which was demonstrated in [7] where the properties of lightest 
mesons in hot and dense environment making use the virial expansion were estimated in chiral 
perturbation theory. In the present paper similar questions are discussed within the instanton liquid 
model of the QCD vacuum [^. 

The quark generating functional for stochastic ensemble of gluon conflgurations being deflned as 
(the sign of double averaging means the corresponding average over the gluon flelds) 

- J Dij^D^ (( e^(^'^''^) ))a , 

at nonzero temperature and flnite chemical potential is estimated by the saddle point method 
and could be presented (up to the details inessential for our consideration here) with an auxiliary 
integration over saddle point parameter A (we are dealing here with SU{3) gauge group and quarks 
of two flavours A^/ = 2) [9] in the following form 




Here ip'^ = {ipf,ipj), f = 1,2 are the quark flelds composed by the spinors of flxed chirality 

ipy = P± ipf, P± = — 2~^' ^ instanton liquid density, p is the mean size of instanton liquid 

pseudoparticles, V4 = TL^ denotes the 4- volume of system and fii, = (0, //). Working in the Euclidean 
space we are summing up over I by introducing the Matsubara frequencies k4{l) = 2tiT (j. + 2) with 
the convention Ak^ = T. Summing up and integrating over spatial momentum components in Eq. ([T]) 
will further be marked as Then the Lagrangian Cint of 4-fermion quark interaction is given with 
the chiral components as 

and for the right hand flelds one should change L ^ R. 

Currently, the corresponding effective Lagrangians for the L n R components are set in terms of 
the quark zero modes. At zero values of temperature and chemical potential the zero modes approach 
[T0| . being combined with the mean fleld approximation [11], is providing us with very reasonable 
phenomenological results. But at flnite temperature the superposition of (anti-)calorons which are 
the periodic solutions of the Yang-Mills equations [I2] looks much more relevant to be used as 

the saturating conflguration of suitable generating functional. At high temperature, meanwhile, the 
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chromoelectrical part of the gluon fleld is expected to develop the mass term Al in the effective 
Lagrangian (the corresponding term in one-loop approximation has been already found out in [T3]). 
It implies the saturating conflgurations should be mainly rooted in the chromomagnetic sector. (In 
the framework of instanton QCD vacuum model the alternative scenario has also been developed 
with assuming the formation of instanton-anti-instanton molecules [H].) The recent lattice studies 
of correlators of the topological conflgurations corroborated [I5] the chromomagnetic character of 
gauge flelds beyond the critical temperature. 

Unfortunately, even the simplest caloron zero mode (in the momentum space) at zero value of 
quark chemical potential /i may not be presented by the special functions and, hence, is impractical 
for the analytical exploration. It is the major reason for us to develop here the approximate procedure 
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in which we substitute the expHcit expression of zero mode ip[A{T, iJ,);T, fj] for its simplified form 
ip[A{0, ii);0, ii] and imply the corresponding corrections could be calculated perturbatively. Such an 
oversimplification results in the L^^^/j^ kernels which are defined by the zero modes (the solutions of 
the Dirac equations with a chemical potential /i) in the form of /ij-functions as regards 

h4{k4, k] jj,) = ^{{k - n - ik4)[{2k4 + in)fi + i{k - n - ikA)f2] + 
+ {k + fx + ik4)[{2k4 + in) ft -i{k + n + ik4)ft]} , 



hi{k4,k;n) 



TT kj 

4e 



2k- n){k- n- ik4)fi + {2k + ijL){k + /i + ik4)f{ + 



+ 
+ 



2{k - n){k- n- iki) - -^(/U + ik4)[kl + {k - fif] 
2{k + ii){k + 11 + ik4) + ^(At + ik4)[kl + {k + jif] 



/2- + 



with i 



, 3 and k — \k\, if one considers the spatial components of 4- vector k^, and 



± 
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h(z^)K,(z^) 



with Jj, Ki {i = 0, 1) as the modified Bessel functions. Introducing in the same time the scalar function 

k- 

h{k4, k; ij) which is related to the three-dimensional components as hi{k4^ k; jj) = h{k4, k; /i) -j^ ,i — 
1, 2, 3 (when it does not mislead we omit the arguments of hi functions at all) we have then 

with S{p; fi) = {p + i^)~ h^{p; /i) , S~^{p; —fi) =h^{p; ~fj){p + ifi)^ , moreover, for the conjugated 

function it is valid hfj, (p; —A'-) = hfj,{p; fi), and e is entirely antisymmetric matrix €±2 — —621 = 1. Here 
and similar symbols are used for the 4-vectors spanned on the matrices where = (i^T, 1) 
and T is the 3- vector of Pauli matrices, p^ = p^T^ and [/ is a matrix of rotations in the colour space. 
Surely, we have the similar relations for the right hand components 



(^f ijf) 



R,8 



{if) 



with the kernel as 



R^J{p,q;pi)^T'{p;pi)e'^' [/j'^e<^"r+.(g; , 



where T{p;ii) = (p + ifj,)^ h~{p;ii) , T~^{p;—fi) =h~^{p', —fJ^ip + ifi)' ■ As far as the vector- 
function h{p) is spanned on the vector p only the components of matrices {p + i/i)^ and h^{p; /i) are 
permutable. As a result it is easy to understand the validity of the following identities 

T{p; //) = S^{p; , T+(p; = S{p; fi) . 

In what follows we omit the ^i-dependence of matrices 5, T, S'"'",T+ because the chemical potential 
enters the matrices 5+, T always with the positive sign only and the matrices 5, T"*" with the negative 
one. 

After averaging over the colour orientation the leading term of the expansion of four-fermion 
interaction contribution looks like the following 



pLiVf) (V'f ^2^2^))^ = (#(P1) 5(P1) 5+(gi) (?l)) (Vr(p2) 5(P2) 5+(?2) ^^2^(52)) - 

- (^^(Pi) S{p,) S+{q2) (92)) (4^(P2) S{p2) S+{q,) (gO) , 



(and, of course, a similar expression for the right hand chiral component). Eventually, the interaction 
term may be rewritten as 

Ant = y {^fiPi) S{p^) S-'iq,) {rt)fr i^riqi)) {i^Wv^) S{p,) S^q^) (r,+),,, - 

Performing the auxiliary integration over the bosonic fields La h Ra the four-fermion interaction can 
be transformed to the Gauss-type integral (analogously for right hand chiral fields and scalar field 

Ra) 

\ (^tL rt (^t^ SS^ r+ i,^) -A (^tL rt - \Ll , 

which is quite convenient to obtain the effective Lagrangian in the terms of hadronic degrees of 
freedom 

L = {l + a + 7])UV, R= {1 + a -r]) V U+ , 

jj- tt^t" (T^r" 

Then the corresponding Lagrangian density incorporating the interaction of quarks and hadrons is 

= - zX i,\p) [S{p) S^{q) e^^P-")^ (I + a{x) + r^ix)) U{x) V{x) P+- 

(3) 

- T{p) T+{q) e*(P-'')^ {I + a{x) - 7]{x)) V{x) U+{x) P_J ij{q) , 

and calculating two well-known diagrams with two external hadron lines one may extract the 
necessary hadronic correlators according to the corresponding effective Lagrangian [10] 



AN ^ '^^ M2(A;) 

vf (c[M^ [(fci + i/i)(A:2 + + M1M2] M{ki, k2)M{k2, h) 
1^ ( \ = -4 A V st ^dkidk2) [(fci + l^l){k2 + - M^M^] M{k,, k2)M{k2, fci 

- -4 , A^T \b idkldk2) [{ki+ip){k2 + ip)+MiM2]M{ki,k2)M{k2,ki 



(4) 



2 J 



^ {dkidk2) [(fci + tfi){k2 + ip) - M1M2] M(fci, fc2)M(A;2, fci) 
'4- (27r)4 [(A;,+^^)2 + Mf][(A;2 + 2/x)2 + M|] 

The momentum integration here includes the 5-function {dkidk2) = dkidk2 S{ki — k2—p) carrying the 
external momentum p. The vertex function M{ki, ^2) = XS (ki) {k2) = Xv{ki, ^2) has the following 
structure 

v{k^, k2) = A + ^B^+^C^+^D ^ ^ , (5) 
\ki\ \k2\ \ki\ \k2\ 

with the functions 

[kM 



A = [kih + {ki + ijj)h4\ \k2g + {k2 + ip)gi\ + \{ki + ip)h - kihA \{k2 + i/i)^' - k2g^ , , , , , , 

\k2\ 

B = [{ki + in)h - kihi\[k2g + {k2 + ip)gA\ , 
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C = -[kih + {ki + i/i)/i4][(^2 + if^)g - hg^l , 

D = [{ki + ifi)h - kih4][{k2 + ifi)g - k2gA\, 

where h = h{ki; jj)^ = h^^iki, fj,), g = h{k2] fi), g^ = hii{k2] fi), and in the shortened form one should 
imply Ml = M{ki, fci), M2 = M{k2, k2). The product of matrices v{ki, k2)v{k2, ki) is spanned on the 
unit matrix because 

ri 1 \ A ^2 T fci X fc2 T 

v{k2, ki) = A - iB — — - - iC — — - - iD 1,11,1 , 
|fei| \k2\ \ki\ \k2\ 

and denotion M{ki, k2)M{k2, ki) means the corresponding scalar to be singled out. Undoubtedly, 
all expressions of Eq.([l|) are approximate and going to improve them one remembers the matrices 
S{ki)S~^{k2) (and T{ki)T^ {k2) as well) should be linked with the corresponding quark propagators 
in between (the vertices might be interchanged only when the momentum equality ki = k2 is valid) 
which finally leads to the additional terms of the following form 

Tr{A + zATS) I75} Ch + + iM^){A - zATS) j^s} (k + ifi + ^^2) , 

where N = B y^+C* -|^^+-D |^"^ |^|fc ^| " ^^^^^^^^ numerical analysis performed shows the additional 

terms proportional to the corresponding traces of matrices of the spin type Sj = ^Eijk 'Jj'Jk are 
negligible. 

Now going to analyse the thermal and finte evolution of quark mass and quark condensate 
we would like to emphasize that we do not expect our approach to be valid in the chiral limit but 
reasonably reliable in the precritical region. As usual the saddle point is determined by the solution 
of the following equation 

,^ dk M^ik) 

At zero temperature and particular value of the chemical potential (/ic — 350 MeV) there appears 
the singularity in the denominator of Eq.([6|) which is rooted in the solution of equation 

fc^ -/i^ + M(|fc|) = 

and quarks are going to fill the Fermi sphere in. If we do not calculate the imaginary part of the 
correlators ^ the corresponding integrals might be calculated as the principal values only. Obviously 
the accurate calculation foresees an analysis of all poles and cuts in the imaginary p4-plane. However, 
for the verteces generated by instantons this task, unfortunately, is not feasible today because of the 
fairly complicated analytical structure of the zero modes in the complex plane (see, for example, 
|16| where unrealistic (low) thresholds of chiral symmetry restoration have been received for the zero 
mode forms considered). 

The masses of a- and vr-mesons are extracted from hadron correlator expansions at small values 
of momenta Q R-k{p) = P-n-p'^ + . . ., Raip) = Oo- + f^ap^ + ■ ■ ■ and = ^ but the vr-meson 
mass is given by the Gell-Mann-Oakes-Renner relation. In our paper the coefficients a and P are 
numerically calculated and the precision which we are able to provide for the meson characteristics 
is at the level of 10% only (it is exracted from the result dependence on the spacing of numerical 
differentiation). It looks quite satisfactory for the qualitative appraisals which here we intend to. The 
typical thermal behaviours of M„ and are depicted in Fig.l for the several values of chemical 
potentials /i = 100, 160, 240, 260 MeV (a polarization was taken trivial |p| = while calculating) and 
the upper curve there corresponds to /i = 100 MeV. The mass of a-meson diminishes because of a 
decreasing and /3 increasing at the same time. The distinctive feature of our approximate calculation 
is an absence of characteristic structure having usually observed at the crossing of the a- and vr- 
meson curves while chiral symmetry is restored |17| . Two curves of the dynamical mass behaviour 
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Phc. 1: The masses of a- (upper solid lines) and tt- mesons (lower solid lines) as the temperature 
functions for four values of the chemical potential = 100, 160, 240, 260 MeV, the upper line 
corresponds to /i = 100 MeV. Two dashed lines show the behaviour of dynamical quark mass 
corresponding to two extreme values of chemical potential (upper line corresponds to larger value. 

M(fc4(0); /Li) for ^ = 100 MeV (lower dashed curve) and = 260 MeV (upper dashed curve) are shown 
on the same plot demonstrating a slow decrease with the temperature increasing. It is interesting to 
notice the cr-meson mass is in the domain of expansion validity with // and T increasing which 
allows us to hope that the results of analysis with an effective Lagrangian proposed might be taken 
not only as qualitative ones. It is well seen that with /i increasing the curves describing the cr-meson 
mass become steeper and steeper and around — 350 MeV where the process of filling in the Fermi 
sphere starts cr-meson is simply degenerated. Moreover the hadronic correlator Raip) = o:a+PaP^ + - ■ ■ 
develops almost constant behaviour at (3a- — > 0. With /i further increasing the correlation function 
of TTa-meson becomes degenerate but its chiral partner cTa-meson appears to be the 'physical' one. 
Fig. 2 shows the correlation functions of Raa{p) ^-nd Rnip) mesons as the functions of of momentum 
|p| (MsB) for three chemical potential values = 0, 510, 680 MeV. 

The curve where M^. = 2Mt^ is depicted (dashed line) in Fig. 3 together with the line = 
(solid curve). Unfortunately, more or less reliable theoretical predictions for their locations on the {/i 
- T')-plane do not exist nowadays. For example, in the Nambu-Jona-Lasinio model these curves have 
the intersection point with the fi = 0-axis around rather 'high' critical temperature and their changes 
while the perturbation expansion in the saturating configuration is applied is of special interest. In 
the model developed, the curves should move down on the plot and if cr-meson appears beyond the 
Ma — 2M^ curve the decay channel cr — > 27r becomes impossible which favours its experimental 
identification. In general we could conclude that our results as to behaviour of cr- and tt- mesons 
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Phc. 2: Correlation functions Raaip) (upper solid lines) and i?7r(p)-mesons (lower solid lines)as the 
momentum |p| (MeV) functions for three magnitudes of chemical potentials /i = 0, 510, 680 MeV. 

on the /i T-plane are in reasonable agreement with the results obtained in other papers (see, 

for example, [18j) although one distinction is obvious. It concerns the quark mass behaviour (see, 
the respective dashed lines in Fig.l). Aiming to clarify this point we have studied the generating 
functional at zero temperature in the mean field approximation as it was formulated in [TT]. Its 
considerable advantage comes from the possibility to calculate the quark condensate (dynamical 
quark mass) without appealing the saddle point equation. However, in this exercise the zero mode 
approximation leads to the unphysical increase of quark condensate when the chemical potential value 
approaches dynamical quark mass magnitude. It seems this result dictates a necessity to take into 
account the terms neglected even though partly (remember, we have made use the approximation for 
zero modes) and to investigate the role of non-zero modes. But we believe studying the mechanism 
of filling the Fermi sphere in gives more illumination to the origin of chiral symmetry restoration in 
this case. 

In our previous paper [H] we were analysing the possible mechanism of mixing cr-meson and 
the phonon-like excitations of instanton liquid and demonstrated that notwithstanding the seeming 
smallness of coupling constant (small change of dynamical quark mass) the effects of mixing could be 
powerful enough. The specific feature of these effects while depending on n and T is the alternation 
of the heavy and light component roles. Initially a-meson is the heavier component but gn heating 
(or compressing) this field loses its mass quickly, as was demonstrated, and at reaching a certain 
threshold magnitude the phonon-like excitations play the role of heavier components. 

We are not preoccupied with the claim for the accurate quantative estimates of particular 
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Phc. 3: The curve of = 2M.„ (dashed line) and the curve of M„ = M^r (solid line). 

masses or widths because we realize the approximate character of perturbative scheme proposed. 
Nevertheless, we believe this approach could be useful to clarify the possibility of continuing the 
estimates obtained to the region of large chemical potentials (larger than dynamical quark mass) 
and the limits of applicability of zero mode approximation. 
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